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Abstract 

Based on the nonequilibrium Green's function technique, a unified theory is developed that 
covers quantum transport and quantum diffusion in bulk semiconductors on the same footing. This 
approach, which is applicable to transport via extended and localized states, extends previous semi- 
phenomenological studies and puts them on a firm microscopic basis. The approach is sufficiently 
general and applies not only to well-studied quantum transport problems, but also to models, in 
which the Hamiltonian does not commute with the dipole operator. It is shown that, even for the 
unified treatment of quantum transport and quantum diffusion in homogeneous systems, all quasi- 
momenta of the carrier distribution function are present and fulfill their specific function. Particular 
emphasis is put on the double-time nature of quantum kinetics. To demonstrate the existence of 
robust macroscopic transport effects that have a true double-time character, a phononless steady- 
state current is identified that appears only beyond the generalized Kadanoff-Baym ansatz. 

PACS numbers: 05.60.Gg,72.10.Bg,72.20.Ht 



I. INTRODUCTION 



Quantum transport in semiconductors has attracted a great deal of interest over the last 
few decades. The transport theory of carriers in bulk semiconductors as developed on the 
basis of semiclassical Boltzmann or balance equations as well as the nonequilibrium Green's 
function (GF) technique is well documented in many review articles and textbooks (cf., for 
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instance, Refs. [I|l2l3[) and has reached a high level of sophistication. Compared with this 
achievement, there are only a few studies on quantum diffusion in semiconductors. However, 
this fact does not imply that the knowledge about quantum diffusion has not contributed 
to the overall picture of quantum kinetics. It is rather the unified description of ballistic 
transport, diffusion, and hopping that provides valuable insight into the general structure 
of the theory. This paper aims at the construction of a general transport theory by fo- 
cusing on common features that emerge from quantum transport and quantum diffusion as 
well as from transport via extended and localized states. A comparative analysis of quan- 
tum transport and quantum diffusion, carried out on the basis of a semi-phenomenological 
approach,-!^ clearly stressed distinct features of the general theory that already appeared 
in early studies of the current density.- The problem is best explained by treating the mo- 
mentum representation of the one-particle transport approach. It is usually assumed that 
due to translational invariance the carrier transport in homogeneous systems is properly 
accounted for by only one A; vector. Consequently, the carrier transport should be governed 
by the distribution function /^(A;|)f:), which is the solution of a quantum-kinetic equation. 
However, completely different results were obtained from a unified approach to quantum 
transport and quantum diffusion on the basis of the conditional transition probability that 
satisfies a quantum-mechanical Bethe-Salpeter equation.-i^ The main quantity in this ap- 
proach is the distribution function /^(fc, K\t) that depends on both wave vectors, although 
the underlying model is completely homogeneous and not strongly affected by special initial 
conditions. This finding is all the more surprising since the vector k, refers to a deviation 
from homogeneity. In fact, it is not the full k, dependence that is needed in the calculation 
of transport coefficients of homogeneous bulk semiconductors. What enters the approach 
is the quantity Vn/^(fc, K|t)|K=o, which is interpreted to be a virtual disturbance of the 
homogeneous system that allows to probe the drift-diffusion response. The surprising result 
that both wave vectors k and k, appear also in the general transport theory of homogeneous 
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systems is confirmed by the unified picture of transport via extended and localized states. 
In principle, hopping and band transport can be described on the same footing so that each 
approach is derivable from the other one in a straightforward manner. This equivalence is 
of particular importance, when transport in biased superlattices is treated. By tuning the 
electric field applied perpendicular to the layers of the superlattice, the character of trans- 
port can be driven from hopping to the Ohmic regime and vice versa. From an application 
point of view, however, it is natural to expect that both approaches are not equally well 
adapted for numerical purposes. 

To provide further arguments for the statement that also the vector k, enters the general 
transport theory, let us focus on the current density j (t) , which is given by the time derivative 
of the dipole operator j{t) = {l/V)dD{t)/dt (with V being the volume of the system). In 
the momentum representation, an equivalent form applicable to the steady state is obtained 



where Veff{k) denotes an effective velocity. In the majority of transport studies, Veff{k) 
is simply given by the drift velocity v{k) = Vk£{k)/h that refers to extended states with 
the kinetic energy e{k). This result applies, whenever the interaction Hamiltonian Hint 
commutes with the dipole operator, which, however, does not always happen. For instance, 
in the theory of small polarons, there is an extra current contribution j ~ [D, Hint]^, which 
even dominates the carrier transport via the effective drift velocity Vefjik). The current 
density of this kind is expressed by the k derivative of the collision integral calculated at k = 
0. Another example refers to the proper definition of the spin current in semiconductors with 
spin-orbit interaction that recently became a subject of a lively controversial discussion.- In 
addition, the exact switching from band to hopping transport and vice versa requires the 
consideration of both momenta k and k in the scattering terms and GFs.- 

The former semi-phenomenological studies provided a unified approach to quantum trans- 
port and quantum diffusion that covers transport via localized and extended states.-"^ It 
remains the task to put this suggestive transport picture on a firm microscopic founda- 
tion. The most attractive starting point for this purpose provides the nonequilibrium GF 
technique that was initiated by Schwinger- and Keldysb^ and that has been employed by 
several authors^^'ii'i^i'i^ii^'i^ to treat quantum transport. With respect to the time depen- 
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dence, such a fundamental approach is able to cope with the two-time nature of quantum 
evolution. The double-time character manifests itself in coupled evolution equations for the 
spectral function and the statistical propagator. Moreover, memory effects appear due to 
the integration over the full time history. Unfortunately, most applications of the nonequi- 
librium GF method disregarded the entangled time dependence by relying on a sufficient 
homogeneity in time, which is exploited by first-order gradient expansions, when 'center of 
mass' coordinates in space and time slowly vary. This approximation has a serious disadvan- 
tage, namely it rules out the full quantum-mechanical character of transport phenomena. 
In contrast to these approaches, we account for the two-time dependence of GFs in an ex- 
act manner and generalize our former unified approach to quantum transport and quantum 
diffusion.' The basic quantity will be the two-time distribution function f^{k,K,\T,t) that 
satisfies a quantum-kinetic equation. An example, given in the last Section, demonstrates 
the existence of a macroscopic phononless quantum transport mechanism that emerges only, 
when the double-time dependence of GFs is properly taken into account. 

II. BASIC APPROACH 

A. Symmetry of the Green functions 

We are going to focus on the basic physics of quantum transport that is revealed by 
controlled approximations in a one-particle picture that disregards the Coulomb interaction 
between carriers. Consequently, the electron-phonon interaction takes over the indispensable 
role of an inelastic scattering mechanism, which is needed when treating the nonlinear high- 
field transport. In spite of this restriction, we will, nevertheless, consider a rather general 
model that is compatible with many studies of quantum transport in semiconductors. 

The main issue in deriving basic quantum-kinetic equations is the full account of sym- 
metries. Although the translational invariance with respect to the spatial and temporal 
ti coordinates is broken, when external electric E{t) and magnetic B fields are applied to 
the sample, there remains a symmetry of GFs that is very important for the description 
of field-dependent nonequilibrium quantum transport.— This symmetry expresses the fact 
that a translation of spatial coordinates can be compensated by the vector potentials A{r) 
and A{t) of the external magnetic and time-dependent electric field, respectively. For the 
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expectation values of the one-particle propagators G^, we have 

with the abbreviations A(t2,ti) = ^(^2) — ^(^^i) S'lid dA{t)/dt = eE{t)/h, where the 
vector potential of the magnetic field A{r) is given in the symmetric gauge. The symmetry 
expressed by Eq. ([2]) favors the utilization of the so-called Wigner transformed GFs denoted 
by that are invariant under spatial translations 

^^(^1,^11^2,^2) = G^(rl,^l|r2,^2)e^"'("^)^^-^"'(*^'*^)(^^+^^)/^ (3) 

with 

G^(ri + r, ti|r2 + r, ^2) = G^Vi, h\r2, ^2). (4) 

or a constant electric field, G^ agree with gauge- invariant GFs (cf., for instance, Ref. 
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18|). Equation ([3]) leads to simplifications that are most effectively exploited in Fourier 



space. Changing the coordinates according to 

R= ^ , r = r2-ri, T= — - — , t = t2 - h, (5) 

we perform a Fourier transformation with respect to the spatial coordinates r and R to 
obtain 

G^k,K\T,t) = J rfVe^*("-^(^'*))"G^(fc + |-iA(r,t)-A(r)|T,t), (6) 

with 

A{T,t) = A(T-^-^-A(T+iy (7) 

The reduction of degrees of freedom and the separation of the momentum k become more 
transparent in the absence of a magnetic field, when Eq. is converted to the form 

G^{k, k|T, t) = 6{k - A{T, t))G^k\T, t). (8) 

This equation enables the formulation of the transport theory on the basis of GFs G^ 
that respect the symmetry in the presence of electromagnetic fields. In our former quantum 
approach to carrier transport,—!^ we profited from this transformation and from the general 
symmetry relation 

G^{k,K\T,ty = -G^{k,-K\T,-t). (9) 

Note that the wave vector k, refers to a possible deviation from homogeneity and is absent 
in the conventional treatment of transport in homogeneous systems. 



B. Dyson equation 



A nonequilibrium system is completely characterized by only two independent two-point 
functions. From a physical point of view, the decomposition of the full matrix of GFs into 
statistical and spectral components is most attractive due to its clear physical interpretation 
and due to the unambiguous separation of their dynamical role. Loosely speaking, this choice 
of GFs makes clear which states are available and how often they are occupied. Having this 
natural decomposition in mind, we derive kinetic equations for the GFs and G^, which 
effectively yield the density of states and the double-time density matrix, respectively. In 
the momentum representation, the coupled Dyson equations have the forn>2i 

ih^ - e{k) + ieE{t)Vk G^{kt \ k't') 
t 

= ±hj dkiS^j dt^E^ikt I fciti)G^(fciti I k't') 
t' 

t' 

+ I dtjl^kt I kiti)G^{kiti I k't') 



(10) 



t 

j dti^^ikt I kiti)G^ikiti I 



All scattering contributions are included in the self-energies E^, while the time-dependent 
electric field is treated in the vector potential gauge. The derivation of kinetic equations 
for the GFs proceeds by well established steps: (i) new wave vectors are introduced by 
the replacement k ^ k + k/2, k' ^ k — k/2 and (ii) Dyson equations written down for 
G^{k, k|T, t) and G^{k, k|T, —t)* are subtracted from each other. This procedure leads to 
exact quantum-kinetic equations. For simplicity, let us shorten the cumbersome calculation 
by focusing on a non-degenerate electron gas, for which the Boltzmann statistics applies un- 
der idealized conditions. As a consequence of this simplification, there is a strong imbalance 
between and G^ in the sense that G'^ is much 'smaller' than G'^ as the statistical prop- 
agator G^ is proportional to the carrier density. Adopting this approximation, we arrive at 
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the following non-Markovian integral equation for the correlation function 

G<(fc,K|T, t) 



+ eik --)- e{k + -)- th—^Vk + th—^V, 
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+G<(fc + ^, ,^ - ,^i|T - t - ti)S>(fc - K,\T + h) j. (11) 

To proceed further, one has to specify expressions for the self energies, which are dictated by 
the underlying model (e.g., small polarons), the treated scattering diagrams (e.g., T matrix 
approximation), and the possible inclusion of the initial conditions.— As an example, let us 
select the self-consistent Born approximation with coupling functions Uf- that are local in 
time 

(fc, k\T, t) = Y^ Uf{k, ki, k\T, t)G^ki, k|T, t). (12) 

fci,A 

For the widespread Frohlich-type electron-phonon interaction, this equation reduces to 

uf{k,k',K\T,t) = Df{k'-k\t). (13) 

The kinetic Eq. f fTTl) together with Eq. f|T2|) are not in the form that is commonly used for 
the calculation of the current density and the diffusion coefficient. Further transformation 
steps have to be carried out without any approximation. An obvious procedure would be 
the full exploitation of the symmetry based on the transformation in Eq. ([6]), which would 
considerably simplify the approach to nonequilibrium quantum transport. A consequent 
exploitation of the Wigner-transformed GFs G^, as applied in our former approach,— i^Sii^ 
leads, however, to a complete disappearance of the wave vector k that plays a fundamental 
role in the unified quantum transport description of both carrier drift and diffusion. There- 
fore, we suggest another procedure, which is more general and which discriminates between 
the states of the system (G^) and their dynamical evolution (G^). The clear distinction 
between the role of spectral and dynamical GFs, which is a generic feature of nonequilibrium 
field theory, is accounted for in the new approach by maintaining the GF G^ in Eq. (fTTj) 
in its original form and by replacing only the GF G^ by the symmetry-adapted partner 



G> defined in Eq. ([8]). Loosely speaking, the remaining k dependence in is needed to 
simulate a virtual probe of the system that reveals its dynamical response. 

The reformulation of the kinetic Eq. ( ITTl) is facilitated by introducing new functions 



and U\ 



G<{k,ii\T,t) = a {k,K- A{T,t)\T,t), 
Ufik, fci, k\T, t) = U}{k, fci, K - A(r, t) |T, t), 



(14) 
(15) 



that account for the internal order of the kinetic equations with respect to k. By a further 
transformation, the kinetic energy of carriers is separated out by replacing the basic GFs 
and G^ through new ones Rr" and R^, which are defined by 

-4/2 / I 



G>{k\T,t) = -iR>{k\T,t) exp 



h 



t/2 



dre k + A{T,t) - 



A{T+i)+A{T-i: 



G {k,K\T,t) = iR<{k,K\T,t) 

t/2 



(16) 
(17) 



exp^ -Jdr 

-t/2 



+£(fc + | + A(T + r)-l 



e(fc-| + A(T+r)-l 



A{T+^-) + AiT-i] 



A{T+i) + A{T~^-] 



The reformulation of the kinetic equation, which is based on Eqs. f fT6l) and f|T71) . isolates 
rapidly varying phase factors and leads to the conventional collision integral that appears 
in transport theory. The final exact reconstruction of the original kinetic Eq. ( ITT]) is carried 
out by introducing the physical GFs and respective coupling terms via 



i?^(fc,K|T,t) = /^(fc ± -A(T,t),K|T,t), for t ^ 0, 



(18) 



u}{k,k',K\T,t) =Uf{k±^A{T,t),k' ±^A{T,t),K\T,t), for t^O. (19) 

The execution of all transformation steps is straightforward and leads to a final kinetic 
equation, in which a time-dependent renormalization of the bare kinetic energy appears on 
the left-hand side. As these single-particle corrections disappear in the limit t — 0, which is 
relevant for the calculation of all transport coefficients, we neglect these contributions. This 
approximation does not significantly affect the general character of our approach. Inserting 
Eqs. ( HM to ( fT9l) into the kinetic Eq. ( ITTi) . we obtain our main general result, namely a 
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quantum-kinetic equation for the carrier distribution function that characterizes the 
statistical properties of the nonequihbrium system 



EE /* 

q,X fci ^-(^ 



>n 1 n\rr , ^ 



dti 



1^ ^ ^ 



x/<(fci, Ac|T + |, t - - g + ||T - i^,t,)Qt,{ki, q, k\T, t) 

The quantities that appear in this kinetic equation are given by 

= A{T+i- t,) - A{T - i), AA_ = A(T - ^ + t,) - A(T - 

Pt, (fci, g, K|r, t) = Ft{k, k\T - ^-)Ft, (fci, k|T + ^ - h)Ft,-t{ki, K\T+'^--h] 
Qi,(fci,g,A€|T,t) = Fi(fc,K|T-^)F4(fci,g,-K|T-^ + ti)F_!(,^_,)(fci,K|T-^+ti), (23) 



(21) 
(22) 



Fj,(fc,g,K|t2) =exp J ^ j dre (^fc - g - | + A(t2 + r) - ^(ta)) >, 



(24) 



ti 



Ft,{k,K\t2) =exp|-^y rfr [£(fc-| + A(t2 + r)-A(t2)) (25) 



+£ (^fc + ^ + A(t2 + r) - A(t2) 



The general result in Eq. fl2U]) has a number of intriguing features. First of all, we mention 
its two-time character. All quantities that enter the kinetic Eq. ( I20l) depend on two time 
variables that are responsible for the evolution on different time scales. This generic double- 
time nature of the nonequihbrium dynamics manifests itself in memory effects,— which are 
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revealed by the non-Markovian time dependence in the colhsion integraL In general, the 
time convolution occurs in both time domains, but it is assumed that the most prominent 
effect primarily happens on the microscopic time scale. For the calculation of transport 
coefficients, only the special one-time distribution function f^{k, K,\T,t = 0) is needed. 
However, its determination from Eq. fl20p is still confronted with the double-time nature of 
the problem that resides in the right-hand-side (RHS) of this equation. A solution of this 
reconstruction problem is offered by the generalized Kadanoff-Baym ansatz.^"^'^^ But in this 
approximation, the two-time dependence is notoriously discarded. The physics that emerges 
beyond this approximation is captured by maintaining the general double-time quantum- 
kinetic Eq. (I2U1). 

Let us add a general remark concerning the interpretation of the double-time kinetics. As 
shown in the next Subsection, the distribution function /^(fc,K|T, t = 0) fully determines 
the temporal order of macroscopic transport phenomena. However, this quantity is only 
given by a strict reference to its two-time extension. Moreover, the self-energy is generally 
constructed from many-particle GFs that depend on several time variables. Therefore, on 
a fundamental level, the change of events cannot be brought into a sequence that can be 
mapped on a one-parameter flow of time extending from the past to the future. In quantum 
statistics, the changeability of things seems to be a more general concept than the temporal 
evolution based on an absolute time in the sense of classical physics. 

Most essential for the construction of a unifled kinetic theory that covers both the carrier 
drift and diffusion is the presence of the second wave vector k in Eq. ( l20l) . which, generally 
speaking, refers to a spatial inhomogeneity. Similar to the role played by the time variable 
t, the full K dependence does not enter the expressions for the transport quantities, but only 
K, gradients at k = (for instance, the basic quantity for quantum diffusion^ is given by the 
vector Vnf^{k, k\T, 0)|k,=o)- Besides the doubling of spatial and time variables, we mention 
an additional peculiarity of the quantum-kinetic Eq. ( !20l) . namely the shift of momentum 
variables in Eqs. (1241) and ( l25l) by the vector potential of the electric fleld. This explicit fleld 
dependence of scattering gives rise to intracoUisional fleld effects and nonlinear transport. 

At the end of this Subsection, we will show that the general two-time quantum-kinetic 
Eq. (!20l) reproduces a number of established results that appear, when further assumptions 
are adopted. The most familiar form of kinetic equations is obtained, when the macroscopic 
time scale is inert to microscopic fluctuations (T + At ^ T in and U^). In this case, we 
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obtain 



oo 

= 5Z / rfti/<(fci,K|T,t-ti)iy(fci,fc,K|T,t,ti), 



/<(fc,A^|r,t) 



(26) 



where the scattering probabihty W is easily read off from Eqs. fl20l) to fl25l) . For a better 
readabihty, the result is presented in the Appendix. This kinetic equation still preserves the 
double-time nature of GFs, but restricts to locality in time with respect to the macroscopic 
time scale T. From Eq. fl26l) . an important sum rule for k = follows. Calculating the sum 
over A;, we immediately obtain from the left hand side of Eq. (126|) the equality 



W{ki, k,K = 0|T, t, ti) = 0, 



(27) 



which is confirmed from Eqs. (l20l) to (l25i) by taking into account the symmetry relation 



U^{k,k',K\T,t) = U^{k',k,-K\T,-t). 



(28) 



The sum rule in Eq. ( 1271) (and its multi-band extension) plays an important role in the 
kinetic transport theory. 

The kinetic description further simplifies, when the double-time character of the problem 
is completely neglected: /^(T, t) f^{T). Specializing to electron-phonon interaction of 
the Frohlich type and restricting to k = 0, we obtain for the scattering probability 



oo 

W{k',k\T) =2Re^ J dtie-'''D>{q\ti)f>{k' + q\T) 



(29) 



X 



[p{k' + f , - f , q\T, h) - p{k' + f , + f , q\T. h] 



with the following field-dependent phase factor 

il 



P{k\k,q\TM) =expi ^y"rfr £(fc' + |+ j dT'F{T')) 



T-tl+T 



T-ti 



T-tl+T 1 V 

-5(fc'-|+ I dr'Fir')) U^^ 

T-ti J ^ 



T-ti 
,fc+ / drFir) 



(30) 



Here, s denotes the Laplace variable of the rudimental microscopic time variation that 
regularizes the ti integral, and F{t) is an abbreviation for eE{r)/h. From Eq. (|29|) . the 
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sum rule in Eq. (127|) is again easily verified. Finally, we obtain the kinetic equation for the 
transport under the influence of a constant electric field (E(t) ^ E) in a form that was 

n 

published many years ago^ (cf., also Ref. (^). 
C. Current density 

An expression for the current density is naturally derived from the general conservation 
law for the particle number. To illustrate the procedure, let us treat the kinetic equation 
for particles with kinetic energy 6{k) that are scattered via the Coulomb interaction. The 
conservation law is easily expressed by nonequilibrium GFs (cf, for instance Ref. js^]), 
and a straightforward calculation leads to the following result expressed in the momentum 
representation 



Scattering does not explicitly enter this equation, in which n denotes the carrier density. It 
is sufficient to calculate the distribution function f'^{k\T,t = 0) = f'^{k,K = 0\T,t = 0), 
which does not depend on the momentum k. How general is this conclusion? Obviously, 
hopping transport is not captured by Eq. fl^ as localized states have no dispersion. In 
fact, Eq. (pri) is a special result applicable to the well-studied models that can be integrated 
into a more general definition of the current density, which covers more complex systems 
with higher-order scattering. According to this definition, the current density is expressed 
by the time derivative of the dipole operator j(t) = {l/V)dD / dt. This physically appealing 
approach is indeed more general than Eq. ( ISTl) . In the momentum representation, we have 



which is converted into another equivalent form by taking into account the kinetic Eq. fl26|) 




(31) 



fc 




(32) 




(33) 



fc 



with an effective velocity given by 



oo 




— oo 



(34) 
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The drift velocity is denoted by v{k) = Vk£{k)/h, and the vector Wi{k,k'\T,ti) is an 
abbreviation for V^W^k, k', K,\T,t = 0,ti)|K,=o. Whenever the sum rule for this vector field 
"^j^Wi^k, k'\T,ti) = is satisfied, Eqs. ( 15^ and (IMl) reproduce the conventional result 
given in Eq. flHTl) . This fortunate situation happens, for instance, for the Frohlich electron- 
phonon coupling and the Coulomb interaction. In general, Eq. fl31l) is applicable, when the 
interaction Hamiltonian commutes with the dipole operator. However, this condition is not 
always fulfilled. For instance, for the transport of small polarons, the Wi contribution in 
Eq. is most essential so that only Eq. (132!) [or the equivalent Eqs. ( |33!) and (134|) ] provides 
meaningful results. The definition of the current density in Eq. ( l32l) includes the k gradient 
of the full distribution function f^{k,K,\T,t = 0) at k = 0. This k, dependence reappears 
in Eq. ([31]) via the vector Wi. Consequently, it is not sufficient to deal with a distribution 
function that depends only on one quasi- momentum k. In fact, the general basis for treating 
carrier transport is provided by the kinetic Eq. (l20l) . from which the k, dependence can be 
determined. 

In summary, we conclude that Eqs. (|33l) and (|34l) put our former semi-phenomenological 
approach^ on a firm microscopic basis and lead to a general expression for the current density 
that takes into account the two-time character of quantum transport. 

To illustrate the additional ability of the approach to simultaneously cover transport via 
localized and extended states, let us, for simplicity, treat the steady-state transport in a 
one-time approximation under the infiuence of an applied electric field Edc- The momentum 
representation in Eqs. ( l33l) and ( IMl) is adapted to the description of low-field transport, when 
the states remain essentially extended. With increasing field strength due to Wannier-Stark 
(WS) localization, negative differential conductivity can appear. To describe this transport 
regime in a more appropriate fashion, the expression for the current density in Eqs. fl33|) and 
is rewritten in an exact manner.- Up to intracollisional field effects, the result for the 
current density along the direction of the electric field 



is compatible with negative differential conductivity j ~ 1 / E^c- Moreover, Eq. fl5Sl) proves 
that there is no current in the absence of any inelastic scattering. A strong electric field 
generates Bloch oscillations that localize carriers so that there is no current as long as only 
elastic scattering is present. By switching to the Houston representation, one arrives at 




(35) 
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another equivalent expression^ 

oo 

J = enYl E (^^)/^(fclX(fcl' fc^)' (36) 

which clearly reveals the hopping character of the transport, d denotes the periodicity of 
the lattice {Id is the hopping length), and A;^ is the momentum perpendicular to the field. 
The /-sum extends over the whole WS ladder. The potential of the approach is illustrated 
by its ability to unify the band and hopping picture by an exact reformulation that mediates 
between them. 



D. Diffusion coefficient 



What favors our approach to quantum diffusion is its close relationship to the carrier 
drift treated in the previous Subsection. Let us follow the same line of reasoning by first 
focusing on the regular part of the diffusion coefficient^ defined by 

Doit) = dV^^{^p\x)^P{x)), (37) 

with X = {r,t) (the spin variable is not indicated). The equivalent expression in the mo- 
mentum representation has a form that is similar to Eq. (J32ll 



k 



(38) 



k;=0 



Again, we mention that both wave vectors k and k appear in this definition. Therefore, 
the basic kinetic equation that describes quantum diffusion has to be formulated for the GF 
/^(fc, K,\T,t), which comprises not only A; but also k. Introducing the vector field 

g{k\T,t) = tf<{k\T,t) = tVJ<ik,K\T,t)U=o, (39) 
we obtain the equivalent form 

{oo 
v{k)-g{k\T,0) + z I dUg{k\T,-h) -J^Wiik^k'lT^h) 
-oo 

oo 

-^E / dtJ<{k\T,-h)J2W2ik,k'\T,h), (40) 



-oo 
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in which W2 denotes the second derivative V^VI^|k,=o- The exphcit scattering contributions 
in Eq. fHOl) indicated by Wi and W2 vanish for widespread models such as the Frohhch 
electron-phonon couphng and the Coulomb interaction. In contrast to the current density, 
which is governed by the distribution function diffusion phenomena are described by 
means of the vector g, which satisfies its own kinetic equation that is easily obtained from 
Eqs. ( 126!) and (139!) . This procedure unambiguously determines the contribution Do(T), in 
which, however, the irregular part is still missing. To complete the calculation of the total 
diffusion coefficient, the vector g in Eq. (HOj) is replaced by a new quantity ip that solves the 
same kinetic equation as g, but with a modified inhomogeneity that is compatible with the 
constraint 

J2v>ik\T,t) = 0. (41) 
fc 

Accordingly, the basic quantity ip, which determines quantum diffusion via the diffusion 
coefficient D{T), satisfies the quantum-kinetic equation 



d dAiT,t) 
dT dt 

+v{k)f<{k\T,t) - 5^i;(fci)/<(fci|r,t) 



00 

cp{k\T,t) = Y^ / dti(p{ki\T,t-ti)W{ki,k\T,t,ti 

1^ J 



(42) 



00 

+2^y" dtif<{ki\T,t-ti] 



Wi{ku k\T, t, ti) - J2 Wi{ki, fcsiT, t, h] 



which is clearly in line with the sum rule in Eq. (HTll . 

To familiarize oneself with the derivation of basic results concerning quantum diffusion, 
let us work with the Carson-Heaviside transformation of the kinetic equation with respect to 
the time variable T [/(s) = s dT exp{—sT)f(T)] in the treatment of a constant electric 
field. From the kinetic equation for g and the sum rule in Eq. fj^Tl) . we obtain 

(43) 



gik\s,t) = cpik\s,t) + -B{s,t)f<ik\s,t), 

s 



in which the following quantities appear 

B{s,t) = Y,'"effik\s,t)f<ik\s,t), 



(44) 



V. 



00 

:k\s,t)=vik)+z I dti:Ci^i|i^5^1¥i(fc,fc%ti)- (45) 
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Replacing the vector g in Eq. (HOj) by the new vector field ip according to Eq. (H3|) . we obtain 



ioo 



oo 

-^E / dtJ<ik\s,-t,)J2W2ik,k'\0,t,) 

-oo f'' 

+-{Y.''eff{k\s,0)f<{k\s,0)\ . (46) 
^ k ^ 

The last term on the RHS of this equation is nothing but the irregular contribution, which is 
subtracted out according to the proper definition of the diffusion coefficient.^ Consequently, 
only the first three terms on the RHS of Eq. (1461) survive and define the total diffusion 
coefficient D{s). These results provide a rigorous theory of quantum diffusion, in which the 
double-time character is accounted for by a vector field that is the solution of the specific 
quantum-kinetic Eq. (l42l) . At this stage, the theories of quantum transport and quantum 
diffusion have reached the same level of sophistication. 

The above theory of quantum diffusion is formulated in the momentum representation, 
which is adapted to extended states. By exact manipulations, other equations for the dif- 
fusion coefficient are obtained that are more appropriate in the WS regime, when carriers 
execute Bloch oscillations. Within perturbation theory with respect to scattering, for which 
the Hamiltonian commutes with the dipole operator, we obtain for the steady state 

D.. = Yl v.ikMk) = E (^(^^) - ^(^^))' f^ik')W{k', fc), (47) 

which has the same structure as Eq. ( l35l) derived in the previous Subsection. In the targeted 
regime of field-induced localization, there is no diffusion without inelastic scattering. 

A more general expression for the diffusion coefficient applicable to the WS regime is 
derived within the outlined approach by exploiting the WS representation.^ The final result 

oo 

^ = 2 5Z E ildff<{k'^)W^;l{k'^, k^), (48) 

fcx.fcY oo 

with W being an effective scattering probability, allows an interpretation within the hopping 
picture that relates carrier diffusion to the square of the hopping length [Id)"^, the lateral 
carrier distribution function f^{k±), and the scattering probability in the site representation. 
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III. EXAMPLE: PHONONLESS TRANSPORT 



Based on the nonequilibrium GF technique, a unified approach has been developed that 
covers both quantum transport and quantum diffusion and that is hkewise apphcable to 
transport via extended and locahzed states. A sahent feature of this theory is the double- 
time character of quantum transport. A natural question arises: What is the significance of 
this double-time dependence? An answer is gained only beyond the generalized Kadanoff- 
Baym ansatz. Generally speaking, it is difficult to draw an overall conclusion regarding the 
physical potential of the two-time quantum kinetics. Summarizing the bulk of conventional 
transport studies, it is tempting to assume that the double-time approach resolves only 
minor corrections that are more or less unimportant. That this assessment cannot be the 
whole truth will be illustrated by a macroscopic transport phenomenon that has no analogy 
in the conventional approach because of its strict double-time character. To be more specific, 
a steady-state current will be identified in the WS regime that is driven by dc and ac electric 
fields without the participation of any inelastic scattering. This phononless current appears 
only beyond the Kadanoff-Baym ansatz and is due to the double-time dependence of the 
GFs. 

A. Solution of the kinetic equation 

The double-time dependence is studied by a model calculation that is simple enough to 
allow for an analytical solution. Some results obtained by a K-independent approach have 
already been published previously.— The model refers to a one-dimensional semiconductor 
superlattice, which is biased by dc and ac electric fields 

E{t) =Edc + EacCos{uja), (49) 

that are sufficiently strong so that WS localization occurs {^Idc^ ^ 1 and ^lac'^ ^ 1? with 
T being an effective scattering time and flac,dc = eEac,dcd/h). Bloch oscillations that appear 
in this transport regime are accounted for by a discrete Fourier transformation of the GFs 

oo 

fHk\T,t)= J2 fHT,t)e''''. (50) 

(=— oo 

To calculate the Fourier components of the double-time distribution function f^(T,t), we 
treat scattering on polar optical phonons with energy hwo and neglect the smooth q depen- 
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dence of the coupling [D<{q\t) D^{t)]. Within the WS regime, only the / = component 
of the GFs enter the collision integral so that Eq. (120|) takes the form 

oo 

l-^ + Un,, + Un^,cos{ujacT)cos('^^ dt, (51) 

k,q Q 

X { [D<{t - ti)<l',,,(fc, q) - D>{Um,^{k, q)] f<{T - |, t - h)f>{T + t,) 

- [D>{-h)<!>l,{k, q) - D<{t + h)^t.,o{k, q)] f^{T - |, t + U)f>{T - -t^)}, 

with the following field-dependent phase factor 

t 



t 

e[k + q + A{T + T-^--U)-A{T-^-)] (52) 



-e[k + A{T + T~\-t,)-A{T-^-, 



(53) 



The double-time character of the approach is still present in Eq. (!5T]) . The non-Markovian 
behavior extends both over the macroscopic (T) and microscopic {t) time scale. The main 
source of the T dependence is the ac electric field that appears directly on the left hand side 
of Eq. (l5Ti) . In most approaches, the double-time dependence is neglected by omitting the 
t dependence in f^{T,t), which results from microscopic scattering processes described by 
the RHS of Eq. (ISTI) . However, both different lines of time evolution are generally coupled to 
each other by a convolution integral, the field-dependent kernel of which determines the role 
played by the kinetic history. An analytic solution of Eq. (!5T!) is found for weakly coupled 
superlattices with the dispersion relation 

£(k) = ^{l-cos{kd)). (54) 

Considering the periodicity with respect to the T dependence 

oo 

ff{T + 2iT/uj,,,t) = ff{T,t), ff{T,t)= J2 ffMe'""^^^^, (55) 

m=—oo 

the main Fourier component f^{m = 0,t) of the WS regime {Qac ^ 1) is calculated from 
the homogeneous integral equation 

oo 

j dt^[D<{t,)f<{0,t^)f>iO,t-t^)-D>{t^)f>iO,t,)f<{0,t-t^)] = 0. (56) 
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This equation, derived under the condition of narrow minibands A/hQac,dc ^ 1, deter- 
mines the time dependence of the distribution function that appears beyond the generahzed 
Kadanoff-Baym ansatz and that was ignored in most previous approaches. A solution is 
searched for in Fourier space by adopting the ansatz 

f<{0,uj) = f>iO,u;)fiu;), (57) 

and by considering the normahzation condition 

oo 

/ ^f^{0,u;)f{u) = l. (58) 

— oo 

Inserting the expressions 

= ihlE^ [^^^ + + - ^^^^ 

for the electron-phonon couphng, it is easily verified that an exponential function in u solves 
the Fourier-transformed version of Eq. fl56l) . To determine the prefactor from Eq. fl58|) . the 
function /q'^(0,u;) is needed. In our previous studies,—'^ we obtained for weakly coupled 
superlattices (A — 0) 

^ f ^V4f/ - 00^, \uj\<2VU , ^ 

/>(0,c.)= <^ ^-^'^ , (60) 

I 0, otherwise 

with f/ denoting the coupling strength of white-noise elastic scattering on impurities. The 
most remarkable features of this density of states function are its nonanalytic character with 
respect to the coupling U and the absence of tails at the band edges. From Eqs. (158|) and 
fl60|l . we obtain the final result for the distribution function 



which applies whenever carriers thermalize more quickly in a given quantum well than they 
need to escape by tunneling. l\ denotes the modified Bessel function. With increasing 
miniband width A, the layers are more strongly coupled to each other, and the solution in 
Eq. (!6T|) is no longer adequate. In this case, only a numerical solution of Eq. (!56l) is available 
that accounts for a non-Markovian time evolution. 
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B. Current density and diffusion coefficient 



Despite previous results concerning the hopping transport in the WS regime, let us look 
for a phononless transport mechanism by exploiting the more general two-time approach. 
The amazing result will be that there is in fact a phononless transport, when the double-time 
dependence of the GFs is properly accounted for. 

Let us first focus on the constant steady-state current that is driven by external dc and 
ac electric fields under the exclusive infiuence of short-range elastic scattering on impurities. 
The steady-state current density 

27r/Wac 

1 =-5: / iTf<(mt = 0) (62) 



1 de{k) ijja 
h dk 27C 

k 

Ad 1 
'2h2i 







en—- [/<_i(m = 0,t = 0) - /<i(m = 0,t = 



is expressed by the components /;^_|_]^(0, 0) of the distribution function that according to 
Eq. (1201) obey the kinetic equation 

1^ + z/fidc + «/^^acCOs(cU,eT)|/<(T,0) (63) 

oo 

= UJ2 / dh [<f ot. {k, q) - {k, q)] 

X {/o<(T - |, -ti)/o>(T - |, ti) - f<{T - |, t,)f>{T - |, Pi{T). 

It is a consequence of this equation that no current can fiow through the superlattice, when 
the variations on the microscopic time scale disappear: /o^(T, t) — > f^{T). This fact confirms 
the general conclusion mentioned in previous Sections that within the one-time picture only 
inelastic scattering enables carrier transport in the WS regime. However, the double-time 
nature of the kinetic evolution opens up a new channel, which enables phononless transport 
of Bloch oscillating carriers. This specific transport mechanism appears only beyond the 
generalized Kadanoff-Baym ansatz. 

The formal solution of Eq. (!63l) has the form 



< ^ ^ Jk^^(l^)jk(l^) 

lYYi — rv~i 1^. — 1-^-. ^ ' ^ ' 



(64) 

— \ LO / \ bJ / I.\Lj„ -I- K.ljJ„„ 

fc= 

where Pi^m. are the Fourier components of Pi{T), the calculation of which follows the same 
steps as outlined in the previous Section. A straightforward procedure applicable to weakly 
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FIG. 1: The normalized drift velocity Vd calculated from Eq. (I65p as a function of ^d^T for = 3, 
iOacT = 2, QacT = Ij and f7r^ = 0.05. The quantity Vdo is given by A^rd/(357r^^), and r denotes 
the scattering time. 

coupled superlattices (A —* 0) leads to the final result for the current density in steady state 

7rUA^dJ,{u')J_,{u') ^ J|(z/') 



J = envd, Vd 

oo 



%h? uo"^^ sin(7rz/) ^-^ k + u 



(65) 



27r 



/o>(0, u;)/o>(0, + cu,e(fc + z/)) [/(cu + cu,e(A; + z/)) - /(a;)] , 



with the abbreviations u = Qdc/^^ac and u' = Qac/^ac- This result again confirms that 
there is no phononless current, when the trivial solution f{uj) = 1 is accepted, which is 
suggested by the Kadanoff-Baym ansatz within a strict one-time approach. The specific 
constant current contribution in Eq. fl65|) disappears also, when the ac field is switched off 
{Qac = 0) because the combined density of states vanishes for 2^ldc > VU- Consequently, 
the assertion that there is no current in the hopping regime without any inelastic scattering 
remains valid also in the two-time approach, when only a constant electric field is applied. 

The calculated drift velocity is shown in Fig. 1 as a function of the dc electric field. 
Photon replicas centered around Qdc = kujac are separated by gaps that result from the 
edge structure of the combined density of states. As a drawback of the simplified treatment 
of scattering, singularities appear in the current-voltage characteristics at Vide = kujac- A 
damping of these resonances is easily accounted for in more refined and realistic approaches, 
which unlikely alter the qualitative physics discussed in this Section. 

The diffusion coefficient is obtained by a similar calculation. According to Eq. (j46l) . D is 



21 



expressed by the function Lp 

1 de{k) uOa 



^-^ilrt I W.«^0,, (66) 

k Q 

which satisfies the following quantum-kinetic equation 

^ + ilVLdc + il^ac cos{uJacT) COS (^^^^ |v5/(^, t) (67) 
{fiUT,t) - fUT,t) - {f<{T,t) - /_<i(T,t))} 



d 



Ah 

2n/d 



J dke-'^^'^Y^ j dtiip{ki\T,t-ti)W{ki,k\T,t,ti) = Pi{T). 



^ 



This equation differs from Eq. (!63|) for fj^ (T, t) by the appearance of a inhomogeneous term. 
The formal solution has the form of Eq. ( !64l) . Restricting the calculation of Pi^m{t = 0) to its 
lowest-order contribution, we obtain the final result for the diffusion coefficient applicable 
to the WS regime 

^ ^ ^ Jj iS^acl^ac) Vdd 

k=—oo 

which can be formally interpreted in terms of a relationship suggested for the current density 
more than 40 years agc^ 



Divide, ^lac) = Jl[ — ] D{^ldc + kuJac). (69) 

Accordingly, the diffusion coefficient under the combined infiuence of ac and dc electric fields 
is easily obtained from the quantity 

/^(l^.e) = ^COthf^V (70) 



2 V A . 

which refers to the absence of the ac electric field Viae = 0. However, this interpretation has 
only a formal character and heavily depends on the approximations made in the derivation. 
The main conclusion is the same as for the current density, namely that this kind of high-field 
quantum diffusion appears only beyond the Kadanoff-Baym ansatz by a strict treatment of 
the double-time dependence. The experimental demonstration of this exclusive double-time 
quantum effect should be feasible by studying biased quantum-box superlattices. 
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IV. SUMMARY 



Starting from a semi-phenomenological kinetic approach, a unified one-electron theory of 
quantum transport and quantum diffusion was developed in previous worka^i^ that applies 
both to the band picture applicable to extended states at low electric fields and to the 
hopping picture for transport under quantizing electric field. Both approaches are completely 
equivalent and can be mutually derived from each other. For the current density, this 
equivalence was already demonstrated in Ref . [2^ . Furthermore, a comparative treatment of 
carrier drift and diffusion^"^ revealed the particular nature of quantum diffusion. Whereas the 
drift velocity goes back to the nonequilibrium distribution function, the diffusion coefficient 
turns out to be constructed from a derived quantity that does not solve the Boltzmann 
equation or its quantum-kinetic extension. Most disturbing was the necessity to deal with 
the total GF /^(fc,K|T, t) that depends on two wave vectors A; and k. The k dependence 
seems to be dispensable for the description of transport in homogeneous systems, which are 
translational invariant. However, to probe carrier diffusion, at least an initial inhomogeneity 
of the carrier ensemble is necessary so that the k, dependence must be preserved in the unified 
description of drift and diffusion. 

In this paper, we put the former semi-phenomenological approach on a firm microscopic 
basis by applying nonequilibrium GF techniques. The unified theory of quantum transport 
and quantum diffusion has been constructed from the quantum-kinetic equation for the 
full distribution function /^(fc,K|T, t). The most salient feature of this extension are the 
appearance of the double-time nature of quantum kinetics and the related non-Markovian 
evolution in two time channels, namely the microscopic and macroscopic time regime. On 
the fundamental microscopic level, quantum evolution seems to be more general than the 
classical schema that dictates a strictly one- dimensional progression from the past to the 
future. An interesting question concerns the relevance of the two-time quantum kinetics, 
namely whether it is possible that new physics appears in this domain. A preliminary 
answer provides the treatment of a one- dimensional superlattice subject to dc and ac electric 
fields. The existence of phononless carrier transport and diffusion is demonstrated, the 
origin of which is the two-time dependence of the GFs. This distinct steady-state transport 
mechanism appears only beyond the generalized Kadanoff-Baym ansatz. Its experimental 
verification seems to be feasible by studying quantum-box superlattices. 
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The rigorous two-time quantum-kinetic approach presented in this paper is hkewise apph- 
cable to quantum transport and quantum diffusion and covers both transport via extended 
states and hopping between locahzed carriers. 
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APPENDIX A: SCATTERING PROBABILITY IN EQ. (f26|) 



Assuming strict locality in time T, we obtain from the kinetic Eq. (l20l) and the definitions 
in Eqs. to 



X 



W{ki,k,K\T,t,ti) 

^|e(ti)/>(fci - g - ^ IT, t,)Pt, (fci, q, k\T, t) 

\_5k^,k+q+AiT+t/2-ti)-A(T~t/2)Ux{ki — q, ki, K,\T,t — ti 



- Ski,k+A{T+t/2~ti)~A{T-t/2)l^x {ki - — , fci - qr - — , 0|T, ti) 



e(-ti)/>(fci -q + ^\T, ti)Qt, (fci, q, k|t, t) 



X 



K 



Skuk+A{T-t/2+h)-A{T-t/2)l^x{kl + —,ki - q + —, 0\T,ti 



— Ski^k+q+A{T-t/2+ti)-AiT-t/2)l^xiki — q,ki, K,\T,t — ti)] |. 

This equation completes the two-time quantum kinetic Eq. (l26!l . 
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